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On the QCD coupling behavior in the infrared region 
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The summary of nonperturbative results for the QCD invariant coupling a s obtained by lattice 
simulations for functional integral and by solution of approximate Dyson-Schwinger equations 
reveals a puzzling variety of IR behavior of a s (Q 2 ) even on a qualitative level. This, in turn, rises 
a question of correspondence between the results obtained so far by different groups. 

We analyze this issue in terms of mass-dependent coupling constant transformations and con- 
clude that the question of the IR behavior of effective QCD coupling and of propagators is not a 
well-defined one and needs to be more specified. 
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1 Introduction 



This paper is devoted to the issue of QCD invariant (effective) coupling function a behavior 
in the infrared (IR) region. 

The notion of invariant coupling is the central one in current practice of quantum-field 
applications of the renormalization group (RG) method. In QFT, the very existence of RG 
is generically connected with Dyson's finite renormalization transformations. The invariant 
coupling (IC) — or invariant charge — has been introduced (see pioneer papers |], ||, |[], 
as well as [|J]) via a product of finite Dyson's renormalization factors Zj , the product that 
is invariant with respect to the mentioned transformations, without any reference to weak 
coupling or ultraviolet (UV) limits. 

By expressing Z{ via values of renormalized QFT functions one usually obtains an ex- 
pression for IC a(Q 2 ) as a function of one argument Q = \/ Q 2 = Q 2 — Ql — momentum 
transfer (or reference momentum) — in terms of a product of the coupling constant, and the 
vertex and propagator Lorentz-invariant amplitudes taken in the momentum representation^] 
in some particular renormalization scheme. 

Such definitions have been introduced and then used in the perturbation case for the UV 
and IR asymptotics in the mid-50s. In current practice, ct(Q 2 ) is usually employed in the 
UV limit only. The massless minimal subtraction renormalization scheme MS turns out to 
be the most convenient one and, hence, more popular in the UV analysis. 

Meanwhile, the mentioned definitions are valid in a more general, mass dependent, case 
- see refs.[||, [| f|, [|. In the gauge QFT models containing several vertices with the same 
coupling constant this mass dependence, in turn, yields a specific vertex dependence, ie. ; the 
dependence on a particular dressed vertex chosen for the IC defining. In QCD, for instance, 
one can use three-gluon, gluon-ghost, four-gluon and various gluon-quark vertices. In the 
massless MS scheme, all ICs constructed with the use of different vertices are the same. 
However, in the massive case, they could differ essentially || in the Q < m region. 

Here, the comment on the difference between massless and mass-dependent subtraction 
schemes and related RG solutions is in order. For example, there exist two variants of the 
MS scheme. In the common one, for the QFT model with one coupling constant g , with 
the help of massless counter-terms one constructs a beta-function /3(g) , the function of one 
argument. The RG equation solving for IC results in expression g = f(£, g) depending on the 
logarithm i = ln(Q 2 / fi 2 ) and describing the UV asymptotic behavior. Meanwhile, in a more 
general formalism, the RG generators (beta and gamma functions) are defined with the help 
of the MS renormalized, mass-dependent approximate expressions. They are functions of 
two dimensionless arguments Q 2 /m 2 and g . The RG equations solving produces expressions 
with the additional "massive" argument m 2 / fi 2 . The general expression for IC is now of the 
form g(Q 2 / fi 2 , m 2 / fi 2 , g) . The mass dependence is essential in the IR domain. This algorithm 
was used in QEDj2|, |7j for describing the IR singularity of the electron propagator. 

This mass dependence, in our opinion, reveals itself in the current analysis of the IR 
properties in QCD, in particular, at lattice calculation of the functional integral and in 
solving the truncated Dyson-Schwinger equations (DSE). Various groups of researchers use 
different definitions of coupling constants and IC in the IR region. For instance, in the 

^ee, e.g., , eqs.(|), (g) and (§) below. 
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IC defining the Tubingen group usesf| the gluon-ghost vertex, the Paris group || [Tcfl 



three-gluon one, while the "Transoceanic team" [11, O - - the gluon-quark vertex 



Section 2 contains a short overview of different groups' results for a s obtained by solution 
of truncated DSE as well as by lattice simulations in the IR region. This summary reveals 
a rather wide variety of IR behaviors. 

In Section 3, we remind the origin of the IC a s (Q 2 ) notion and note that it is free of 
any relation with the weak coupling or UV massless limits. Then, we consider the item of 
nonuniqueness of the a s (Q 2 ) definition, and the related question of coupling constant and 
IC transformations. We discuss this issue in terms of coupling transformations similar to 
the renormalization-scheme ones in the weak coupling limit, but, generally, more involved 
and mass dependent. Here, the important feature is the vertex dependence that is typical of 
QFT models with gauge symmetry. 

Further on, in Section 4, we return to the item of correspondence between different IR 
behavior of a s (Q 2 ) . In particular, we give few examples of transformations that are singular 
in the IR region and, in some cases, results in the a s (Q 2 ) behavior in a similar way to the 
ones discussed in Section 2. 

Our conclusion is that the question of the IR behavior of the effective QCD coupling 
a s (Q 2 ) (and of propagators) is not a well-defined one. The results of different groups, for- 
mulated in the momentum representation, should not be compared directly with each other. 
Calculation of hadronic characteristics remains as a reasonable criterion for comparison of 
different lattice simulation schemes. 



2 Nonperturbative results in the IR region 

Here, we give a short overview of nonperturbative results obtained by different groups in 
calculating the effective coupling ct s {Q 2 ) behavior in the IR region by numerical simulation 
on lattice and by approximate solution of the truncated Dyson-Schwinger equations (DSEs). 

Each of these groups uses its own definition for a s (Q 2 ) , compatible with the asymptotic 
freedom regime in the UV region, and obtains its own result in the IR domain. 

2.1 Tubingen: gluons and ghosts 

We start with a description of recent results obtained by the Tubingen group. This group 
activity (for a review see paper ||) proceeds in parallel along two lines: solving the truncated 
Dyson-Schwinger equations and the lattice simulation of Euclidean 4-dimensional functional 
integral for quantum gluodynamics. Thus, this group limits itself to the gluon-ghost sector 
of QCD that turns out to be convenient for the DSEs truncation in the Landau gauge. Here, 
the gluon-ghost vertex renormalization function V due to the Ward-Slavnov identities can 
be represented via analogous functions of the gluon Z and ghost G propagators. At the 
same time, in this gauge, V drops out from the invariant coupling defined as 

a Tu (Q 2 )= a s Z(Q 2 )G 2 (Q 2 ). (1) 

Due to this, only two propagators enter into the properly truncated DSEs. 
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Lattice simulations for propagators reveal [13|, [14| a rather specific IR behavior. The gluon 
function Z(Q 2 ) - the transverse amplitude of the gluon propagator — passes through 
maximum at Q = ^/Q 2 ~ 1 GeV and then goes to zero approximately as Q , the ghost one 
G(Q 2 ) is monotonous and has a power IR singularity close to Q^ 1 ^ 2 , while the product ([I]) 
tends to a finite value. This picture is supported by 

DSEs solving that yields the power IR behavior for the propagators 
G{Q 2 ) ~ (QY K ; Z(Q 2 )~(Q 2 Y with re = 0.595 
and finite limiting value for the invariant coupling 

a Tu (0) ~ 2.97(= 8.92/iV c ) . 

2.2 Paris group, 3— gluon vertex 

Here, the invariant QCD coupling function is constructed on the basis of the three-gluon 
vertex with the nonsymmetric MOM subtraction in the Landau gauge for the QCD model 
with two quark flavors. In this particular MOM scheme the invariant coupling 

a P {Q 2 ) = a s t 2 {Q 2 )Z\Q 2 ) ; f(Q 2 ) = T 3gl (Q 2 ,0,Q 2 ) , (2) 

according to lattice simulations, obeys a very peculiar behavior in the IR region. It passes 
through the maximum at Q 1 — 1.5 GeV and then quickly approaches zero at Q ~ 0.5 GeV 
- see Fig. 3 from M. 

The data here are not sufficient to define the IR limiting regime. Nevertheless, in our 
opinion, the power decreasing ~ (Q 2 ) u ; v > 2 is not excluded. 

2.3 "Transoceanic team", gluon— quark vertex 

This teamf] considers (for a recent review see ]12[) the QCD model with two massive quarks. 
It defines its invariant coupling gxo on the basis of the gluon-quark vertex in the particular 
MOM scheme with gluon momentum equal to zero. It reads 

9to{Q 2 ) = a s T TO (Q 2 )Z 1/2 (Q 2 )S(Q 2 ) ; r TO (Q 2 ) = IV„(0; Q 2 , Q 2 ) , (3) 

with S , the quark propagator amplitude. The results for ^xo obtained by lattice simulation, 
qualitatively, are close to the "Paris group" ones - see Fig. 4 in paper JTlJ] or Figs. 8 and 9 
in |T2[ . However, the IR asymptotics of IC, evidently, has a different power 

g(Q) ~ {Q 2 Y' 2 ■ a^{Q 2 Y- M < 1 . 

2 The group includes the authors from various centers from Australia, Great Britain and continental 
Europe. 
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2.4 "ALPHA" — The Schrodinger functional 



The "ALPHA" group considers the QCD model with two massless flavors. It uses the 
Schrodinger functional (SF) defined in the Euclidean space-time manifold in a specific way: 
all three space dimensions are subject to periodic boundary conditions, while the "time" 
one is singled out — the gauge field values on the "upper" and "bottom" edges differ by 
some phase factor exp (—77). Then, the renormalized coupling is defined via the derivative 
r" = dr/drj of the leading part of the effective action 

r = a _1 r + t 1 + a s r 2 + . . . 

as (cf. eq.(2.43) in ref . [ 15 1 and (8.3) in ref. ||16||) a function in the coordinate representation 

a SF (L) = ryr' (4) 

with L being the spatial size of the above-mentioned manifold. 



To follow the «sf(£) evolution, a special trick was used fig] . The "step scaling function" 

a(a SF (L)) = a S F(2L) 

has been introduced. For explicit implementation of a , the beta-function is necessary. 
On the other hand, numerically, this function can be defined from lattice simulations by 
comparing the results for lattice L with 2L — see Fig. 16 in ref.[|TH]. Meanwhile, fresh results 



17|j reveal the steep rise of the SF running coupling in the region asF — 1 and L — > oo 



Here, the analytic fit for the numerically calculated behavior of «sf has an exponential form 

a SF (L)~e mL with m~2.3/L max (5) 

and L max , a reference point in the region of sufficiently weak coupling, indirectly defined via 
the condition asF(-^max) = 0.275 . For a physical discussion of the IR (i.e., at Q 2 — > ) QCD 
behavior, in the papers of ALPHA group the "usual quantum-mechanical correspondence" 
L — 1/Q is used. 



3 Mass— dependent coupling function 

Renormalization schemes with scale parameter /i coinciding with subtraction momentum 
- the so-called MOM schemes — are singled out from a formal point of view. For these 
schemes the condition of the invariant coupling normalization is of a simple form 

g\ Q 2 =fl 2 = g„. (6) 

In a more general caseQ, that is relevant to minimal subtraction schemes, in particular to 
the MS one, this condition contains the "normalization function" 

9\ Q ^ = N(g) , (7) 

3 In this section we shall use a general notation g for a coupling constant along with a more specific one 
a = g 2 /Air commonly adopted in gauge theories. 
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that can be dependent on mass(es). 

The normalization (|6|) corresponds to the simplest functional equation (FE) for the in- 
variant coupling, that in the massless case is of the form 

_ /x _ \ Q 2 
g{x;g») = g [^;g(t;g^J , x = — . (8) 

Here, the condition (H) follows from the RG functional equation (|8|). We shall refer to 
as to the "canonical" normalization condition. Meanwhile, for the case (0), the FE is of 
a more involved form — see below eq. ([l4]) depending on the mass argument m as well. 



3.1 Massive case 



Now, we consider a more general case that takes into account an important generalization: 
the g dependence of mass(es) of particle(s). Here, an invariant coupling acquires one more 



argument that can be chosen (see, e.g., ref.||]) as a ratio m 2 //i 2 
condition for an invariant coupling 

g~N(i,y;g) = N(y,g) 



y , and the normalization 



(9) 



generally differs from the simplest, "canonical" , form @ . It contains a function N depending 
on two arguments. 

Note here that we have introduced a new special notation g^ for the coupling function 
with the general normalization condition (|9|) to distinguish it from a coupling function nor- 
malized to g . In what follows, we use the term "invariant coupling" (IC) for the function g^ 
with the nontrivial normalization function N ^ g and the term "effective coupling" (EC) for the 
function g with (||). 

The effective massive coupling satisfies Q a rather simple functional 



g(x,y;g) = g ~,g{t,y;g) 



_ Q 2 m 2 

' ~ ~7i ' y ~ ,,2 ' 



(10) 



and differential equations 

dg(x,y;g) 



x- 



dx 



x 



P ( -> 9 J , x p-g(x,y;g) = o 



with 



d d n , . d 
ox ay og 



P{y,g) = 



(12) 



The last two equations (|TI]) and ([12]) are generalizations of massless equations, and turn to 
the latter in the limit y — > . We define a special notation Xp for the particular infinitesimal 
(Lie) group operator (that in the given context was first introduced by Ovsyannikov fT8[| in 



1956). Its superscript refers to the coordinate of the last term. In what follows we shall refer 
to X as to the Lie-Ovsyannikov (LO) operator. 
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The effective coupling g by definition satisfies the canonical normalization condition 

g(i,y;g) = g, (13) 

while other RG one-argument functions, like propagator amplitudes, at Q 2 = y 2 are nor- 
malized to unity. 

At the same time, the functional equation for the IC is of a more involved form 



9N(x,y;g) = gN 



x y 



n 



y 



9N(t,y; g) 



(14) 



that reflects the group composition property and corresponds to the normalization condition 
([]). Here, n is the function reverse to N with respect to the second argument. 
By differentiating (|i~4"D in two various ways one obtains the differential equations 



dg N (x,y;g) 
dlnx 



B ( - ,g N {x,y;g) 
^x 



with 



d _ 

B(y,g) = J7Z9N {£,y;n{y,g)} 



Xp ■ 9N(x,y,g) = 



l=N(y,g) 



(15) 



(16) 



Note here that for IC we have two different beta-functions B and (5 entering into eqs.([15|) 
and related by the second of equations flilf ). 

Meanwhile, it turns out to be possible to treat invariant coupling g^ as an effective 
coupling function corresponding to some other coupling constant 7 = N(y, g) . Indeed, by 
introducing a new coupling function 



V, 7) = gN{x, y; n(y, 7)} 



(17) 



we find that it satisfies the group functional equation of a simple type (|i~0D , the differential 
eqs. analogous to fllTD with the LO operator Xb with the same generator B{y, r )) and can 
be considered as EC corresponding to the coupling constant 7 . 

On the other hand, if we define a new auxiliary function g by the relation 



g{x, y,g)=n{-, g N (x, y; g)\ = n {- , 7 (x,y; N(y, g)) 



x 



X 



(18) 



it will also satisfy the simple functional equation (|TUp and can be treated on the equal footing 
with EC 7 . Both group differential eqs. for g contain the same generator f3 . We shall refer 
to it as to the "group effective coupling" to stress that, generally, it is not related to IC by 
the Dyson transformation. 



3.2 Mass dependent transformations 

Consider now, in the general mass dependent case, the transformation from the coupling 
constant 7 and IC 7 to the new coupling constant g and EC g . The transition from 
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popular massless MS to the above-mentioned massive MS can be represented as a sequence 
of such transformations. Let, without loss of generality, coupling constants g and 7 be 
connected by 7 = N(y, g) , g = n(y, 7) . Then, according to eq . fll8|) , the effective couplings 
will be related by 

7(x,y;i) = N{y/x,g(x,y;g)} (19) 



(20) 



and the corresponding group generators by fllq). 

In other words, it is possible to consider the transition 

g _> 7 = N(y,g) 



as a coupling constant transformation related to the change of the renormalization pre- 
scription. That is, starting with the set {g,g(x,y;g),[3(y,g)} we can move to another one 
{7, j(x, y; 7), B(y, 7)} by the transformations (p0|) and also by fll9"|) represented in the form 



i(x,v,i) = N \ - ,9 



X 



(21) 



that is a composition of fll?D and (0) with the generators /? and S related by Eq. (p~6j) . 

Correspondingly, in the group FEq., for an RG covariant function s(x,y;g) (like, e.g., , 
propagator amplitude) it is necessary not only to take into account relation (|TJ) between IC 
g~N and EC g 

( \ s(t,y;g) /x y ' 

s{x, y; g) = — —- rr -si-,-; g(t, y; g) 



s(l,y,g(t,y;g)) 



t 1 t 



(22) 



but also normalization of s itself 

V, 9) = S(y,g) ^ 1 . 
The related differential equations are of the form 
d\ns(x,y;g 



dlnx 



7 



-,9(x,y,g) ; X p ■ \ns(x,y;g) = 7 ( - , g 
x \ \x 



(23) 



with 



d 



t=i 



d d 

and j(y,g) = j(y,g) + { y-^ - (3(y,g)— [> In S{y,g) . (24) 



In a given gauge QFT model with one coupling, both EC and IC, besides the renormalization 
scheme, are specified by the choice of the vertex and the way of its subtraction. 



3.3 The "vertex" dependence 

In the gauge theories for the mass dependent case in the definition of effective coupling 
there appears a new specific degree of freedom related to the existence of several Lagrangian 
structures with the same coupling constant. Here, in defining IC, different vertices can be 
used. This means that in such a case we have various possibilities for defining IC in the same 
RS. The transition between diverse ICs is described by transformations just considered. 
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To illustrate, take, e.g., a few dressed QCD vertices, the 3-gluons T 3g i , the gluon-ghost 
Fgi- g h and the gluon-quark T g i^ q ones. 

Generally, each of them could be used for the QCD invariant coupling defining, e.g., 

a 3gi = a s Tl gl Z 3 , a,gi- gh = a s V 2 gl _ gh Z G 2 , a gi - q = a s T 2 gl _ q Z S 2 (25) 

with Z, G and S — gluon, ghost and quark propagator scalar amplitudes. 

In the UV massless limit, all definitions are, in a sense, equivalent: they coincide in the 
MS scheme and could be "slightly different" in MOM schemes due to various definitions of 
the "reduced" vertex functions presented as functions of one space-like argument Q 2 . 

However, in the massive case they could be drastically different in the IR region due to 
diverse dependence on the light quark masses. As it is well known, the quark (fermionic) 
propagator as well as the gluon-quark vertex are singular on the quark mass shell. For the 
light quarks this leads to the singular IR behavior. 

The gluon-ghost vertex and related coupling a g i- g h seems to be less sensitive to the 
mass effects. Indeed, at the one-loop level, quark mass effects can contribute to it only via 
the gluon propagator factor Z . Meanwhile, the polarization loop is free of afore-mentioned 
singularity. Two others ICs, a^gi and a g i- q are more sensitive. In each of them, besides 
three propagator amplitudes, the mass effects come also via the vertices T^i or T g i- q . 

Quite probably, just to this there corresponds a more "quiet" variant of the IR behavior 
for a{Q 2 ) obtained by the Tubingen group. 



4 The model coupling transformations 

Here, we present some model transformations of the coupling constant and related IC trans- 
formations, including the mass-dependent ones. It is shown that an appropriate a s trans- 
formation could drastically change the properties of the corresponding running coupling as 
a function of Q 2 in the IR region. 



4.1 Examples from APT of massless transformations 

Let us start with examples induced by the Analytic Perturbation Theory (APT) that has 
been devised recently |19|, |(J to clean out the perturbative QCD (in a "bloodless" way — by 
imposing the Kallen-Lehmann analyticity) of unphysical singularities like the Landau ghost 
pole. 

In the APT, the transition from the usual invariant MS coupling constant a s to the 
Minkowskian % and Euclidean ole ones can be treated [21j as a coupling transformations 



similar to a change of the renormalization scheme (RS). At the one-loop case 

.1 11 
ct s — > &m\ a s ) = —r~ arccos — -j^^^=^^^= = — — arctan(7r/3 ct s ) , (26) 
ttPo a/1 + -n 2 l3l a s 2 ^Po 

a s ^a E (a s )= a s + -j- (l - e 1 ^ " ) _1 , (27) 

Po 
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Here, the first transition "looks quite usual" as q>m can be expanded in powers of a s 
(like in the usual RS transformation), while the second one in the weak coupling case is close 
to the identity transformation as far as the nonperturbative term ~ e~ 1 ^ /3 ° aa leaves no 
"footsteps" in the power expansion. 

If one starts with the common one-loop a s (Q 2 ) ~ 1/ \n(Q 2 / A 2 ) with its first order pole at 
Q 2 = A 2 , then, clS £L result of substitution of ( |26| ) and (^) into the RG differential equation 
and its integration, one arrives at the ghost-free expressions 



arccos 



VL 2 



Q!an(Q 2 



1 

A) 



L>0 
1 

I 



arctan(7r/L) 
A 2 



A 2 



Q 2 - A 2 



(28) 
(29) 



that are well known in the APT. 

The transformation functions «m(«s) and oie(oi s ) obey important properties. They 
tend to a s in the weak coupling limit 



and 

Besides, they are finite 



tti — ► a s at a s C 1 
are finite at a q = oo . 



a,i — > 1/Po m the limit a s — >• — . 



[AF] 
[GhF] 



The first property [AF] provides a correspondence with the weak coupling limit with its 
asymptotic freedom property. The second one [GhF] reflects the absence of ghosts and the 
[IRf] one relates to the finite IR limit. 

This last feature is absent in the transformation 

a s -> a N = a s (l - e~ 1/A)C * s ) (30) 

that is equivalent to the one obtained p3| in the modified APT. 

Just due to the [AF], [GhF] and [IRf] properties, both APT ICs (H) and (||)0 

(x) enjoy the Asymptotic Freedom, 

(xx) are free of ghost singularities at Q 2 = A 2 ; 

(xxx) have a finite IR limit at Q 2 = . 

4 The same function of the initial coupling constant F(g) ~ arctan(g) appeared [Q in the exact solution 
of the two-dimensional Thirring model. 

5 For both a and a an the corresponding beta-functions have zero at a = 1//3q and are symmetric under 
the reflection [a — l/2/?o] — > — [a — 1/2(3q] . Meanwhile, the beta-function for a(s) turns out to be equal to 
the spectral function for a a n(Q 2 ) • This last property, that provides a peculiar realization of the Schwinger 
hypothesis is valid |2l], [2^] outside the one-loop approximation. 



10 



At the same time, the last transformation ([30]) satisfying only [AF] and [Ghf], yields an 
expression obeying the singular IR behaviour 

Q 2 - A 2 

" N(g2) = /3 Q 2 ln(Q7A 2 )' (31) 

It has an "extra Q~ 2 factor" that, as some people believe, relates to the linear growth of the 
interquark potential. 

Consider one more massless transformation 



e^-expjte- 1 / 2 ^} m 
a s -> «sf = ol s - E — - ; b = — , (32) 



leading to the expression 

, L A 

«sf(Q 2 ) = - ; /02x • — h 7- (33) 



with the exponential singularity olsf ~ e M ^ as Q — > . It reminds an analytic approxima- 
tion (|5|) of recent ALPHA group numerical results transposed (see, e.g., Fig. 3 in Ref-JH]]) to 
the IR momentum region with the help of "quantum-mechanical correspondence" L — > 1/Q . 
Such an essential singularity contradicts the Kallen-Lehmann representation. Note also that 
expression (|3"2"D, in turn, obeys an essential singularity at a s — > +0. 



From the given examples, there follows a simple "rule of correspondence" between trans- 
formations (p6"D — (|32| ) and resulting expressions ( p8|) — (|33| ) for invariant couplings: 



under the coupling constant transition a s — > oti = fi(a s ) the invariant coupling 
function transformation is 

a s (Q 2 ) -> «,(Q 2 ) = f\(a s (Q 2 )) (34) 



4.2 Massive transformation 



Consider now the mass-dependent transformation of a coupling constant. Generally, accord- 
ing to Q2"0|), it involves a function of two variables 

(35) 



a s — ► a* = N(y, a s ) . 
Here, the IC transformation looks like 



a s — > a [x, y; a 



x 



x,y;n\-,a 
x 



(36) 



The r.h.s. of the last relation, in distinction to (34), contains the variable y/x = m 2 /Q 2 
that can influence the IR behavior. 
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To illustrate, take a model expression 



N(y,a) = a M {a s ) 



1 + y 



ca M (a 3 ) 



(37) 



that yields 



a*{Q 2 ) = a{Q 2 



Q 2 - m? 



2\ ca(Q 2 ) 



(38) 



Q 2 + m 2 

Here, the coupling constant om( oi s ) and EC a(Q 2 ) are defined by expressions of the type 
(p6|) and (^3), or by their more involved two- or three-loop counterparts — see, e.g., EH], [25[ . 
Here, it is essential that <5(Q 2 ) is free of unphysical singularities and monotonically increases 
to a finite IR limit a(0) ~ 1 . The structure of mass singularity in (|38| ) resembles the well- 
known IR one (see, 0]) of the QED vertex on the mass shell. 

For positive c ~ 1 and small quark mass m < A values the function d*(Q 2 ) , defined 
by can be very close to the numerical results for a s , obtained by the Paris group (see 



subsection 2.2) and by the "Transoceanic team" (subsection 2.3) 



5 Discussion 

As it has been mentioned above in Section 2, the results of various groups for the QCD 
invariant coupling in the momentum representation — obtained by lattice simulations of the 
Euclidean functional integral as well as by approximate solving of truncated SDEs — turn 
out to be quite different in the IR region. At the same time, their "physical" results for 
hadronic properties of matter seem to be more correlated. 

Our "model constructions" of Section 4 demonstrate that the IR properties essentially 
depend on a precise way of the IC defining. We used a class of ICs that in the UV region cor- 
relates with the perturbative QCD a s coupling in the MS scheme. Just this correspondence 
is usually considered to be essential in lattice simulations and analysis of SDEs solutions. 

In practice, these "admissible" ICs, satisfying the condition [AF], can correspond to 
diverse lattice calculations. As far as these calculations satisfactorily describe confinement 
and hadronic physics, it is reasonable to consider them as "physical" ones. 

Special attention should be paid to EC g 2 {L) of the ALPHA group with its exponential 
growth ~ e ML with lattice spatial size L . To interpret this result in the momentum transfer 
representation, one needs to be very cautious with performing the Fourier transformation as 
far as the usual Tauber criterion is not valid. This issue will be considered in more detail 
elsewhere. 

To conclude, we believe that there is no direct physical sense in attempts to establish 
some "correct IR behavior" of the perturbative QCD invariant coupling. Any "infrared QCD 
physics" like hadronic and r decay ones, cannot be described in terms of only pQCD notions. 
Either pQCD description has to be supplemented by some additional semiphenomenological 
parameters like "effective parton masses" p6| , p7| , [28|j (of order of pion mass) and anomalous 
vacuum averages, or some other, intrinsically nonperturbative, means, like lattice simulations 
or the Dyson-Schwinger equations, should be used. 
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